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ABSTRACT 

Aims. We present an analysis of archival Spitzer images and new ground-based and Hubble Space Telescope (HST) near-infrared (IR) 
and optical images of the field of M83 with the goal of identifying rare, dusty, evolved massive stars. 

Methods. We present point source catalogs consisting of 3778 objects from Spitzer Infrared Array Camera (IRAC) Band 1 (3.6 
pm) and Band 2 (4.5 pm), and 975 objects identified in Magellan 6.5m FourStar near-IR J and images. A combined catalog of 
coordinate matched near- and mid-IR point sources yields 221 objects in the field of M83. 

Resuits. We find 49 strong candidates for massive stars which are very promising objects for spectroscopic follow-up. Based on their 
location in a fi - F versus V - I diagram, we expect at least 24, or roughly 50%, to be confirmed as red supergiants. 
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1. Introduction 

Massive stars are of prime importance in astrophysics. The con¬ 
sequences of their birth and evolution have profound ramifica¬ 
tions for galaxy evolution. Throughout their lives, their physical 
characteristics shape and mold their immediate environment via 
intense radiation fields and strong winds. Many massive stars 
expire as supemovae (SNe), thus affecting subsequent star for¬ 
mation and the chemical enrichment of galaxies, their interstellar 
media, and the intergalactic medium. Predicting the evolution of 
a massive star relies heavily upon the mass loss experienced by 
a star over the course of its life. 

Mass loss in massive stars (M > SMq), particularly in the 
late stages of evolution, is poorly understood (see the recent re¬ 
view by Smith|20T4 l. Wind mass loss rates of massive stars on 
the main sequence, once thought to be a major contributor to 
massive star mass loss, are now being revised to be 2 to 3 times 
lower than previous estimates based on observations in the ultra¬ 
violet and optical (see Crowther et al.|20(j2} Repolust et al. 2004 


Puls et al. 2006[ for example) as well as in the X-ray (|Cohen et al. 


201 l| l. Episodic and eruptive mass loss in evolved massive stars 
must therefore be more important than previously thought. 

Understanding mass loss via observations of Galactic mas¬ 
sive stars is difficult. Massive stars are rare and are formed in 
regions along the Galactic plane where extinction can be patchy 
and very high. Eor example, [Ramirez Alegrfa et al. (2012 1 esti¬ 
mate visual extinctions up to 21.5 mag for members of the young 
massive stellar cluster Masgomas-1 at a distance of ~3.5 kpc. In 


* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 

** Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 


order to observe massive stars at large distances in the Milky 
Way, it becomes necessary to observe them in the infrared (IR). 
By far the best efforts to discover new Galactic massive stars ex¬ 
hibiting heavy mass loss have been most successful using Spitzer 
24 pm observations that revealed circumstellar shells of ejected 
material surrounding massive evolved stars (e.g. Gvaramadze 
|et al.|2010{ [Wachter et al.|2010j l. In order to form a more com¬ 
plete census, however, we must look outside our own Galaxy. 

The consequences of late evolutionary stage mass loss for 
lower mass stars may be seen in the progenitors of some super¬ 
novae. One specific example is SN 2012aw, a type IIP event in 
the galaxy M95 ( Fagotti et al.||2012 i. Pre-explosion mass loss 
is likely responsible for the dusty circumstellar material (CSM) 
seen as a significant visual extinction around the red supergiant 
(RSG) progenitor that has a mass at the high end of the range 
for type IIP SN progenitors ([Eraser et al.||2()I21 [Van Dyk et al. 

[2m2| ). 

The situation concerning mass loss in more massive stars 
is uncertain. SN 2009ip was initially mistakenly identified as 
a SN ( Maza et al.||2009| l based on its initial rise in brightness. 
However, it never reached typical SNe magnitudes, and was 
eventually understood to be an object undergoing a luminous 
blue variable (LBV) eruption ( Berger et al.||2009| Miller et al. 


20091. SN 2009ip also experienced a similar outburst in 2010, 


followed by a brighter event in 2012. The jury is still out con- 


|Smith||2014[ Mauerhan et 

al.[[2014t [Graham et al.[[2014[l or an 

imposter ([Pastorello et al. 

2013 

Margutti et al. 

2014 |Fraser| 


et al.|2015| l , but poorly understood mass loss events have clearly 


played a vital role in shaping the late stages of evolution for this 
50 - SOMq star. The specifics of episodic and eruptive mass loss 
mechanisms in both the higher mass (> 3QMq, [Humphreys & 
Davidson 1994| Humphreys et al.|2014[l and l ower mass regimes 

' imith[[2' 


(RSGs with masses < 30Mq, [Smith[|2014|l remain unknown. 
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While the case of SN 2009ip is rare, discovering and cataloging 
all massive stars with mass loss is a tractable investigation. 

Mass loss in massive stars can reveal itself via an IR excess 
in different ways depending on the type and mass of the star. In 
a recent study of luminous and variable stars in M31 and M33, 
[Humphreys et al.| ( |2014[ l found that LB Vs had no hot or warm 
dust evident in the IR. Feii emission stars have an IR excess 
attributed to circumstellar dust and nebulosity, both of which 
are indicative of mass loss. Meanwhile, only some A-F spectral 
type supergiants show IR evidence of mass loss, possibly be¬ 
cause some are post-RSG objects, while others are still evolving 
to the red on the Hertzsprung-Russel (HR) diagram. For redder 
objects (RSGs), mass loss is seen as thermal-IR excess from hot 
dust and molecular emission ( |Smith||2()14| l, and masers in the 
extreme cases Pabing|1996] . 

Work categorizing the massive stellar content of nearby 
galaxies in the IR with Spitzer began with the Large Magellanic 
Cloud (LMC; [Blum et al.||2006l [Bonanos et al.||2009] l. Small 


(Thompson et al. 


Magellanic Cloud (SMC; Bonanos et al. 2010 1 and M33 


2009 1 ). [Bonanos et al. ' p009| [201^ laid the 


foundation for differentiating the types of massive stars by look¬ 
ing at objects with previously known spectral classifications in 
the LMC and SMC. [Thompson et al.[ ( [2009[ ) focused on the 
dust-enshrouded massive stars in M33 identifying several can¬ 
didates as extreme asymptotic giant branch (EAGB) stars that 
exhibited very red colors from Spitzer InfraRed Array Camera 
(IRAC) Band 1 3.6 pm and Band 2 4.5 pm observations, and 
are likely similar to the progenitor of transients like SN 2008S 
( [Prieto et al.|2008| l. [Khan et al.[ ( [2010| ) extended this effort to the 
nearby galaxies NGC 6946, M33, NGC 300, and M81. Of partic¬ 
ular interest was the identification of the brightest mid-IR point 
source in M33, dubbed Object X ( [Khan et al.|201 l[ l, a likely bi¬ 
nary system composed of a massive O star and a dust enshrouded 
red supergiant ( Mikolajewska et al.|2015| l. A similar search was 
also applied to look for rj Car-like objects in NGC 6822, M33, 
NGC 300, NGC 2403, M81 NGC 247 and NGC 7793 ( [Khan] 
et al.|201? I. Britavskiy et al. ([2014|l used Spitzer photometry to 


select bright candidates in the Local Group dirr galaxies IC 1613 
and Sextans A, and their follow-up spectroscopy discovered five 
new RSGs and one yellow hypergiant candidate. 

We aim to continue to characterize the massive stellar con¬ 
tent of nearby galaxies with high star formation rates starting 
with M83. M83 lies outside of the Local Group at a distance 
28.4 mag, [Herrmann et al.[[2008] Radburn 


of ~4.8 Mpc (p __^_ 

Smith et^[2011 1. It is a face-on galaxy classified as SAB(s)c 


([de Vaucouleurs et al.[199l|) with the hfth-highes t Ha luminosity 


in the local 11 Mpc’ volume (log Lh„ =41.25, 


Kennicutt et al. 


2008[l. The star formation rate for a galaxy is directly propor¬ 


tional to its Ha luminosity ( [Murphy et al.|20lT| l. M83 has also 
been host to six historical SNe, five of which have been identi- 
hed as massive star or core-collapse; the Type IIP SN 1923A 
( Lamplarid] 1923[ [Rosa & Richter [1988[ l, the unclassihed SN 


945 B ([Lilier|[ 199011, the TypeTTSN 1950B ( [Haro & ShapT^ 


9501 [Weiler et al.f[T986ll, another T ype II SN1957D ( [Kowal 


& Sargent[[1971||Weiler erar][T986ll, the Type IIP SN 1968L 


Wide Field Camera 3 (WFC3) observations in order to investi¬ 
gate the massive stellar content of M83. 


Kim et al. (20121 studied the stellar content of select regions 


of M83 with a subset of the WFC3 data used here and came 
to the conclusion that young stars are more likely to be found 

( fT999l l 


in concentrated aggregates along spiral arms. Larsen 


reached a similar conclusion, noting that highly crowded clus¬ 
ters are less than 10 Myr old, roughly the expected age for RSGs 
of masses similar to that of the progenitor of SN 2012aw. Also, 
Chandar et al. ( |2010 1 , using HST observations, noted that M83 
contained a large number of clusters with M > 10^ Mq and 
T < 10^ yr. We therefore focus on methods attempting to find 
massive star candidates that are demonstrably separate in the 
data from cluster candidates. We also target the spiral arm re¬ 
gions of M83, in order to investigate the youngest stellar popu¬ 
lations. 

Section 2 describes the observations and data reduction, with 
some initial analysis. Section 3 describes the combined catalog 
between data sets from section 2, and our methods for select¬ 
ing massive stars from the data. In section 4, we compare our 
method of selecting evolved massive stars to other methods of 
photometric selection and present conclusions in section 5. 


2. Observations and Data Reduction 


2.1. Spitzer Observations 


( Bennett|1968[ B^bon et al.|1979[), and the Type Ib SN 1983N 
( Thompson et al.|1983 Porter & Filippenko|1987| l. The high Ha 
luminosity coupled with the copious amount of historical SNe 
and SN remnants ( Blair et al.|2012 2014| l imply that M83 has 
a large population of massive young stars. These stars evolve 
relatively quickly, meaning there may also be a large number 
of dusty evolved massive stars that are prominent in the high- 
quality Spitzer archival images. We also supplement Spitzer data 
with ground-based near-IR and Hubble Space Telescope {HST) 


We extracted the mosaic image of M83 from the Local Volume 
Legacy Survey (LVL, Dale et ak][2009| l in all four of the IRAC 
bands. The hnal image analyzed here covers an area of 15' x 15' 
with a pixel scale of 0'.'75 pix ' centered on the nucleus of M83 
(J2000.0: a =13:37:00.9, 6 = -29:51:56). 

Initial point source lists in Spitzer IRAC 3.6 pm (Band 1) 
and 4.5 pm (Band 2) bands were extracted using the IRARH 
implementations of the DAOPHOT ( [Stetson 1987 \992\ suite 
of programs. A point-spread function (PSF) was built for each 
bandpass from bright, isolated stars. Aperture photometry of 
the PSF stars was converted to the Vega system using the 
aperture corrections and zero-point fluxes given in the IRAC 
Instrument Handboolj^ PSF magnitudes for the entire set of 
point sources were then calibrated by applying the offset be¬ 
tween PSF magnitudes and corrected aperture magnitudes for 
the PSF stars. The coordinates for point sources in the 3.6 pm 
list were checked against those in the 4.5 ;um list using J. D. 
Smith’s IDL code match_2d.pro with a half pixel search 
radius (O'.'375). Positional coordinates were derived from the 
World Coordinate System (WCS) information contained in the 
header of the archival LVL image for the 3.6 pm frame. 

To eliminate the possibility of including foreground or other 
previously identihed objects, we checked the positions of candi¬ 
dates in our list against a number of catalogs. Mimicking the 
method used to match sources between the 3.6 pm and 4.5 
pm bands, we compared our candidate positions to those in 
the PPMXL catalog ( [Roeser et al. 2010p , discarding 73 can¬ 
didates having both a position match, and a proper motion in 
either right ascension or declination of > 15 mas yr * (un¬ 
certainties in proper motion measurements were on the order 
of 10 mas). Larsen ( 1999|l studied young massive star clusters 


' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 


for Research in Astronomy (AURA) Inc., under cooperative agreement 
with the National Science Foundation. 

^ http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandboo 
^ http://tir.astro.utoledo.edu/jdsmith/code/idl.php 
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in M83, publishing a catalog of 149 clusters. Using the same 


half pixel search criteria, we found matches to 5 of the Larsen 
( 1999| l clusters, and removed them from further analysis. We 
also checked the cluster list from a study of the central region 
of M83 by Harris et al. (2001[), but found no matches. Herrmann 


& Ciardul 


(2009|l list 241 photometrically identihed plane¬ 
tary nebulae in M83 from their investigation into the mass-to- 
light ratio of spiral galaxy disks. A coordinate check against 
this list resulted in no matches with our candidates. In a re¬ 
analysis of European Southern Observatory (ESO) Very Large 
Telescope (VLT) imaging data covering one 6!8 x 6'8 held in 
M83, [Bonanos & Stanek ( |2003[ ) discovered 112 Cepheid vari¬ 
ables, none of which matched the coordinates of any of our 
candidates. The supernova remnant population of M83 has been 
studied in some detail. We compared Spitzer candidates against 
the object lists in |Blair et al. P012|l, re jecting 3 objects, and the 
list of new objects in |Blair et ak j2014l, hnding no matches. 

The hnal 3778 objects in M83 are included in Table [T] The 
positions of sources in RA and Dec are given in degrees in the 
hrst two columns of the table. Based on comparison with posi¬ 
tions of 2MASS stars in the observed held, we estimate the as¬ 
trometry to be accurate to about 0'.'3. The remaining columns are 
made of pairs; the source’s photometry in a bandpass followed 
by the uncertainty in the measurement. The columns start at the 
lowest wavelength, with J, and continue with (described in 
|2.2| for the near-IR), 3.6 pm, 4.5 pm, 5.8 pm, and 8.0 pm (de¬ 
scribed in ^3.1[ ). Sources with no detection in a particular band¬ 
pass are left with a blank entry. Some sources were only detected 
in the near-IR or mid-IR and thus only have measurements in J 
and Ks or 3.6 pm and 4.5 pm. 

We estimated the photometric completeness of the list in a 
method similar to |Bibby & Crowtha~ (2012 1 . We placed detected 
sources in 0.25 mag bins and ht a power law to the bright end of 
this distribution. A completeness level limit of 100% is assumed 
to be where this distribution begins to drop from this power law. 
The 100% completeness magnitudes are therefore estimated to 
be [3.6] 18.3 mag and [4.5] ^ 17.5 mag. These limits corre¬ 

spond to absolute magnitudes of Mp 6 ] ~ - 10.1 mag and Mp 5 ] ^ 
-10.9 mag using our assumed distance modulus of p ^ 28.4 
mag. Therefore, we are only able to completely sample the most 
luminous massive stars in these bands. 


2.2. FourStar Observations 

Observations of M83 were made in J and with the EourStar 
instrument ( jPersson et al.||20T3] l attached to the 6.5m Baade 
Magellan Telescope at Las Campanas Observatory, on UT date 
2014 Jan 2. Data reduction was performed with the ESRED 
packag^ courtesy of A. Monson. The resulting 7-band mosaic 
(MJD=56659.31938) was constructed from 7 dithered exposures 
yielding a maximum exposure time of 335 seconds, while the 
corresponding /fs mosaic (MJD=56659.33165) was the combi¬ 
nation of 10 frames giving a maximum exposure time of 582 


both aperture correction and photometric calibration. Because 
the photometric standard stars were on the same frame as the 
candidate objects, the transformation equations simplihed to a 
zero-point for each band and color term between the two bands. 
The computation of transformed magnitudes thus required a 
source to have a measurement in both the 7 and bands within 
a matching radius of 2 pixels corresponding to ~ 0'.'32. The 
2MASS stars in each bandpass were also subjected to the same 
transformation, and the comparison of our measured values and 
those published in the 2MASS catalog are shown in Table]^ The 
average difference in the 7-band is 0.047 + 0.039 mag, while the 
average difference for the (7 - K^) color is 0.056 + 0.045 mag. 
The average difference in 7 is nearly insignificant, however, the 
similar difference in the (7 - color may indicate some sys- 
tematics with the 7-band measurements. 

Artihcial star tests were used to estimate the photometric un¬ 
certainty of a measurement similar to the method employed by 
Grammer & Humphreys (2013|l. Using the PSE constructed for 


each bandpass image, we added one thousand artificial stars (a 
number roughly equal to the total number of point source detec¬ 
tions matched between bands) to each frame via the ADDSTAR 
routine in the DAOPHOT package. We then performed photome¬ 
try on these targets in the same manner used on the point source 
candidate list and repeated this process three times. We com¬ 
puted the root-mean-square (rms) of the differences between in¬ 
put and output magnitudes for 0.2 mag bins. This rms difference 
was then added in quadrature with measurement uncertainties 
from ALLSTAR falling into that particular bin, to make the fi¬ 
nal estimate of measurement uncertainty. It should be noted that 
this rms difference was typically very close in value to the uncer¬ 
tainty estimates from ALLSTAR. The final point source list of 
975 objects are included in Table[T] Sources detected in the near- 
IR are sorted by RA in column 1, with column 2 listing the dec¬ 
lination, both in degrees. Again, from comparison with 2MASS 
sources, we estimate that the EourStar astrometry is accurate to 
roughly (X^. Columns 4 through 7 correspond to 7 magnitude 
and measurement uncertainty, and /Tj magnitude and measure¬ 
ment uncertainty for sources detected in the near-IR. Eor sources 
with no 7 and Kg measurements these columns are left blank. 
We did not correct the point source photometry for foreground 
Galactic extinction. However, for completeness, the current best 
estimate using the Schlafly & Einkbeiner| ( 20lT[ ) infrared-based 
dust map gives 0.047 and 0.020 mag in 7 and respectively. 

we found the 


Using the same method as described in §2.1 


100% completeness levels to be 18.0 mag (My -10.4) in 7 and 
16.6 mag (M;f^ a; - 11 . 8 ) in K,^. 

2.3. HST WFC3 Observations 


seconds. Astrometric calibration employed SWarp (Beitin et al. 
|2002^ software and stars from the 2MASS ( jSkrutskie et al.|2006| l 
catalog. The final images have a field of view of 11' x 11' cen¬ 
tered on the nucleus of M83 with a pixel scale of O'.'16 pix"^ 
Photometry for the 7 and bands was conducted in the 
same fashion as for the Spitzer data. The PSE varied across the 
four CCDs making up the EourStar instrument, so several PSE 
stars in each quadrant of the final mosaic were selected. The 
PSE stars were chosen from the 2MASS catalog and used for 


Blair et al. (2014 1 studied supernova remnants in M83 with 
seven fields of Hubble Space Telescope {HST) Wide Eield 
Camera 3 (WEC3) observations in multiple bandpasses. Two 
fields come from the Early Release Science Program (ID 11360; 
R. O’Connell, PI) with the remaining five coming from the cy¬ 
cle 19 HST General Observer program 12513 (W. Blair, PI). 
These seven fields cover the majority of the bright disk region 
of M83. Eor the reduction and discussions concerning the pho¬ 
tometry of these data, we direct the reader to [Blair et al.|(|2014 


and references therein). We specifically used the imaging anc 
photometry in the E336W (Johnson U), E438W (Johnson B), 
E555W (Johnson V) or E547M (Stromgren y, easily converted 
to Johnson V), and E814W (Johnson I) bands. The quality of 


http://instrumentation.obs.camegiescience.edu/FourStar/SOFTWARlj^8di^jtion.html 


the astrometry for the HS T data are better than O'.' 1 (Blair et al. 
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3. Results 

3.1. Combined Catalog of Point Sources 

In order to characterize the point sources, sources detected in 
both the near-IR and mid-IR were matched by position as de¬ 
scribed above. Final coordinates were adopted from the shortest 
wavelength, J, owing to the band having the better spatial res¬ 
olution compared to longer wavelength bands. To supplement 
these measurements, aperture photometry using the IRAF task 
APPHOT/PHOT was performed on the remaining IRAC bands, 
5.8 fj.m and 8.0 fim at the pixel positions of sources having J, 
Ks, 3.6 lum and 4.5 /urn photometry. The aperture photometry in 
these two bands was calibrated in the same way as described 
for the PSF stars in the 3.6 fim and 4.5 fim bands in the previ¬ 
ous section, with appropriate constants applied from the IRAC 
instrument handbook. 

This near- and mid-IR combined catalog of 221 point sources 
in M83 is presented in Table 

3.2. Candidate Massive Star Collection 

Identifying massive star candidates among our point source 
catalog presents several challenges. For example, the average 
FWHM of objects in the 3.6 fim band is ~ 1"7, which corre¬ 
sponds to 40 pc at the 4.8 Mpc distance to M83. Thus, virtually 
all the clusters identified in [Larsen] ( |1999| l will be unresolved on 
ground-based images, leaving open the possibility of confusing 
massive stars with compact clusters or unresolved associations. 
We also expect contamination from foreground Milky Way halo 
objects like M-dwarfs, and background galaxies. One way to 
overcome this limitation is to compare the location of objects 
on various color-magnitude diagrams (CMDs) and other photo¬ 
metric diagnostic plots. Another way to address this issue is to 
use higher resolution Hubble images, as employed in the current 
paper. 


3.2.1. Mid-IR Selection Criteria 


We applied mid-IR criteria to select massive stars from our sam¬ 
ple following the Spitzer Space Telescope Legacy Survey known 
as “Surveying the Agents of a Galaxy’s Evolution” (SAGE) stud- 
ied both the EMC (jMeixner et al.|2006| and SMC (SAGE-SMC; 
Gordon et al.|2011|l allowing a characterization of the IR prop 


erties of the massive stellar content of both galaxies (jBonanos 
|et al.pOO^ |2010|). In both the LMC ( jBonanos et al.|2009| l and 
SMC ( Bonanos et al.||20T0 l, massive, evolved stars with previ¬ 
ously known spectral types were shown to group in certain re¬ 
gions of the [3.6] versus [3.6]-[4.5] CMD, thanks to the well- 
studied nature of the brightest stars of these nearby galaxies. 
Specifically, looking at Eigure 2 of both [Bonanos et al. ( |2009| l 
and Bonanos et al. ( [2010| l, RSGs are clumped in a region where 
-0.3 < [3.6] -[4.5] < 0.0 and -13.0 < M [3 6 ] < -9, while super¬ 
giant Be (sgB[e]) stars are located in a band further to the red, 
with [3.6]-[4.5] between 0.6 and 0.9 and -13 < Mp.g] < -9. 
RSGs appear “blue” owing to the suppression of flux in [4.5] 
from the CO and CO 2 molecular bandheads ( [Verhoelst et al.| 
[2009| l. LBVs occupy the [3.6]-[4.5] color space between these 
two groups. There are always exceptions to these generaliza¬ 
tions, and one concern may be that the cuts for RSGs may 
not include stars with dusty envelopes. There is evidence that 
these photometric cuts are leaving out RSG. Eor example, in 
Britavskiy et al. ( 2014| l, the RSG IC 1613 1 was originally 
thought to be an LBV based on its color: [3.6]-[4.5]= 0.67. 


Eollow up spectroscopy, however, revealed it to be an early M 
supergiant. 

Eigure 1^ shows the [3.6] versus [3.6]-[4.5] CMD of point 
sources from Table Sources are plotted as black dots except 
for the two regions empirically shown to contain RSGs and 
sgB[e] stars as discussed above. The RSG region is outlined in 
a red box, and point sources are shown with red dots, while the 
sgB[e] region follows the same prescription, but with the color 
blue. Applying these regional criteria for massive star candidates 
yields 638 RSGs and 363 sgB[e] candidate stars. Also plotted in 
Eigure [2] are prominent m assive stars from the literature: rj Car 
(Humphreys & Davidson 1994|l, M33 Variable A ([Humphreys 


[et al.|2006] l, and Object X ( [Khan et al.|201l| l. To represent the 
accuracy of the photometry, we have plotted ellipses on the left 
side showing the typical uncertainties in [3.6] and [3.6]-[4.5] 
for each one magnitude bin. 

One major difference between the LMC and SMC versus 
M83 is that background extragalactic sources become more nu¬ 
merous in the mid-IR CMD at the distance of M83. Ashby et ak] 
( |2009) l showed that for the Spitzer Deep, Wide-Eield Survey 
(SDWES), galaxies make up most of their detections. This is 
especially true past [3.6] >16 mag, where many of our massive 
star candidates exist. We determined the approximate contam¬ 
ination by extracting 42,465 sources from the final version of 
the SDWES catalog ( [Kozlowski et al.[[2010| l from a randomly 
selected one square degree area. These are plotted in Eigure 
as small grey dots. Eigure is the same as Eigure except it 
also shows the possible regions of contamination by background 
sources from the SDWES catalog. This is not entirely represen¬ 
tative of the expected contamination in the field of M83, as it 
contains all objects within a one square degree region while the 
analyzed M83 Spitzer image covers only 15' x 15'. 

To further explore the background contamination, the same 
selection criteria used to find massive star candidates in M83 
were applied to the sample from the SDWES, yielding 4309 ob¬ 
jects in the RSG area and 5348 objects in the sgB[e] region. We 
then scaled these numbers for the area encompassed by the ex¬ 
tracted Spitzer image of M83, 15' x 15' or 0.0625 square de¬ 
grees. The final estimates for contamination are 270 (42%) in the 
RSG region, leaving 368 RSG candidates and 335 (92%) in the 
sgB[e] region, with only 28 objects not likely to be background 
contaminants. [Khan et al.| ( [2013[ l used similar methods to esti¬ 
mate their contamination, finding it to be much less for galax¬ 
ies closer to the Milky Way. In a study of the asymptotic giant 
branch (AGB) and super-AGB (SAGB) population of the LMC, 
[DeU’Agli et al.[ ( 2014| l showed that the most massive (6 - 8Mq) 
AGBs, the OH/IR stars, exist at luminosities above the classic 
limit for AGBs of L ~ 5 x IO'^Lq or M/,oi = -7.1 mag. The AGB 
region outlined in [Khan et al. ( [2010[ l lies in a slightly reddened 
part of the CMD, redward of [3.6]-[4.5]=0.1 mag and where 
^[ 3 . 6 ] < “10 mag. However, given the photometric uncertain¬ 
ties at the distance of M83, there is the possibility that some 
photometrically identified evolved massive star candidates may 
actually be the most massive and most luminous AGB/SAGB 
stars. 

The strongest massive star candidates from the Spitzer pho¬ 
tometry are those objects that lie outside the background con¬ 
taminated region shown in Eigure These include a number 
of objects that may be very red objects [3.6]-[4.5] > 1.5 mag 
or objects with colors bluer than RSGs. The strongest RSG 
candidates lie in the region bounded by -0.3 < [3.6]-[4.5] 
< -0.1 mag for 15.4 < [3.6] < 16.5 mag and in the bounding 
box above the line connecting the two points of [3.6] = 16.5 
mag and [3.6]-[4.5]= -0.1 mag with [3.6] = 18.0 mag and 
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[3.6]-[4.5]= -0.3 mag. It should be noted that photometric un¬ 
certainties may shift some candidates into and other objects out 
of this region. Regardless, these criteria resulted in the identifica¬ 
tion of 118 candidates from the Spitzer photometry. Further re¬ 
finement was made via plotting the positions of these candidates 
over an image of the galaxy. The extent of the images used in 
this analysis (15' x 15') goes well beyond the disk light of M83. 
While massive stars exist in the outer regions of M83, they will 
be more numerous particularly compared to background galax¬ 
ies in the regions of disk light of M83. Because we want to maxi¬ 
mize our potential for selecting massive star candidates for spec¬ 
troscopic follow-up, we therefore removed detected sources ly¬ 
ing outside of a square region centered on the nucleus of M83 
of 800 pixels or 10'. This region is selected based on the extent 
of the disk of M83 as seen in the 3.6 pm image. The remaining 
RSG candidates within this region are 68 point sources selected 
only from Spitzer 3.6 pm and 4.5 pm observations and denoted 
by "'Spitzer RSG Candidate” in column 15 of Table 


3.2.2. Near-IR Selection Criteria 

In order to expand the massive evolved candidate star list, we 
used the FourStar data in /-band along with the Spitzer 3.6 pm 
measurements to construct a plot of [3.6] versus 7-[3.6] similar 
to Figure 3 in [Bonanos et al.| ( |20091 l. Figure shows our ver¬ 
sion of the plot. We have regions of interest outlined in red for 
the location of spectroscopically known red supergiants in the 
LMC ( [Bonanos et al. 12009) 1, 1.0 < 7-[3.6] < 2.0 and fainter than 
= -12, while the blue outline shows the location of known 
sgB[e] stars, 2.0 < 7-[3.6] < 4.5 and fainter than M[ 3 . 6 ] = -12. A 
total of 148 candidates were found in these regions. At the risk 
of removing possible interesting massive star candidates (mas¬ 
sive stars outside of this region due to binary interactions or SN 
kicks), but increasing our chances of selecting individual mas¬ 
sive stars, we again chose to exclude any objects lying outside 
a box of 10' centered on the nucleus of M83. When combined 
with those selected from our analysis of Spitzer observations 31 
sources were selected using both sets of criteria, leaving 185 
massive star candidates. 


3.2.3. High Angular Resolution HST/WFC3 Optical Selection 
Criteria 

Cross referencing coordinates of the 185 candidates selected 
from near- and mid-IR photometric criteria with positions in the 
WFC3 images, we visually inspected the area surrounding each 
candidate and classified it based on its appearance. We ranked 
each candidate based on the probability of conhrming their mas¬ 
sive nature via spectroscopic follow up of the target. Rank 1 
is used for objects that are the most isolated, with few nearby 
stars of comparable brightness. Those objects with slightly more 
crowded fields, but maintaining the appearance of a single star 
are put into rank 2. In rank 3, fields get more crowded, making 
identification of the individual objects more difficult. Some ob¬ 
jects classified as rank 3 show a diffuse nature. For objects in 
rank 4, crowding was worse than in rank 3 or objects are even 
more difficult to recover. Rank 5 is reserved for the worst cases 
of crowding, diffuse nature of the primary candidate, or a com¬ 
plete lack of recovered object at the position. 

The WFC3 fields do not cover all of the disk of M83. 
Because of this only 112 of the 185 massive star candidates may 
be checked by visual inspection. This resulted in the identifica¬ 


tion of 55 (49%) candidates with ranks 1 and 2 for best suited 
for spectroscopic follow-up observations. 

Figure shows examples of each rank, 1 through 5, while 
the online hgure set shows a postage stamp for each of the 112 
candidates grouped by rank. These ranks, as well as a short com¬ 
ment to explain the rank for each candidate, are given in Table 
for all 185 candidates selected from IR data. Objects not in the 
HST held of view have a comment stating as much. Also listed 
are broadband photometric measurements and uncertainties in 
Johnson U, B, V and I as well as the available near- and mid-IR 
photometric measurements. The uncertainties of the HST opti¬ 
cal measurements are merely the output from DAOPHOT and 
have not been more closely checked. They are listed here to give 
the reader a general idea of the accuracy of the HST photome¬ 
try. The magnitudes listed in Table are not corrected for either 
Galactic extinction, nor reddening internal to M83. However, 
for reference. Galactic extinction is Au - 0.288, As - 0.241, 
Av - 0.182, Al - 0.1 mag from the [Schlahy & Finkbeiner 
( 201 l| l IR-based dust map, and the average M83 internal ex¬ 
tinctions from |Kim et al.| ( |2012| ) based on HST observations of 
the central region are Au - 0.696, Ab - 0.559, Ay - 0.441, 
Aj - 0.256 mag. To illustrate the locations of these candidate 
RSGs, Figure shows the greyscale IRAC Band 1 3.6 pm im¬ 
age used in the analysis with the sky positions of the candidates 
overplotted. Candidates of the best rank (ranks 1 and 2) are plot¬ 
ted as red circles, while those of lower quality (ranks 3 through 
5) are represented by blue circles. 


4. Discussion 


There have been several surveys for massive stars that have used 
broadband optical photometry (e.g. Massey||1998[ Levesque & 
Massey 2012[ Grammer & Humphreys |2()13| l as a principal 
means to find red, blue, and yellow supergiants. Thus, a com¬ 
parison of optical versus IR selection techniques is worthy of a 
brief discussion. 

The primary reason to utilize mid-IR observations to ex¬ 
plore the massive stellar content of a galaxy is that evolved mas¬ 
sive stars, with substantial dust or circumstellar material, will be 
more easily detectable in the IR as opposed to the optical. Deep 
optical photometry of galaxies has not been conducted in a sys¬ 
tematic fashion beyond the Local Group Galaxy Survey effort of 
Massey et al. ( |2006| l. Plus, mid-IR selection criteria have been 
shown successful in the Local Group after spectroscopic follow 
up ritavskiy et al.|2014l l. 

One drawback with the current mid-IR surveys is the dif¬ 
ference in resolution compared to near-IR and optical images. 
Massive stars are born, and mostly die, in clusters or associ¬ 
ations, meaning spatial resolution at the distance of M83 is a 
prominent concern. The PSF of the near-IR observations is just 
below 1" corresponding to a size of 23 pc at M83’s distance. 
This size is smaller than the 40 pc size of the Spitzer PSF, but 
may still easily include clusters, associations or several stars. OB 
associations can have half-light radii of a few parsecs and extend 
over several tens of parsecs, while young massive clusters can 
exist in a volume less than one cubic parsec ( [Clark et al.|20051 l. 
One such example is the cluster Stephenson 2, which contains 
roughly 30 RSGs in a radius between 3.2 and 4.2 pc ( jNegueru^ 
et al. 2012| l. The regions of M83 covered by the seven HST fields 
have an order of magnitude better resolution than the FourStar 
observations. This PSF corresponds to ~(y.'l or just over 2 pc. 
While some objects may still be compact clusters, we can elim¬ 
inate many objects that are larger than 2 pc in size, increasing 
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our odds of observing only individual stars with follow-up spec¬ 
troscopy. 

The second issue with using tnid-IR data concerns the lim¬ 
iting magnitude attainable. We found absolute magnitude limits 
of M [3 6] - -10.1 mag for the completeness of the survey, mean¬ 
ing we are completely sensitive to only the brightest RSGs at 
the distance of M83. Similarly, |Khan et ar](|2010 1 noted that for 


Humphreys 


Grammer & 


Grammer & 


To aid in illustrating where conhrmed RSGs reside, in Figure 
l^we have plotted spectroscopically conhrmed RSGs (grey dots) 
from the Local Group that have B, V, and /-band measurements: 
119 from M31 that are not labeled as “crowded” objects ([Massey 


et al. 2009| l, 46 from M33 that are not “multipl e” ([Massey et al. 
2006 1 , 73 from the Milky Way (Ducati 2002[ Levesque et al. 


NGC 6946 at a distance of 5.6 Mpc, the 3cr detection limit in 
the 4.5 fim band corresponded to an absolute magnitude of close 
to -10. Therefore, given the limits in both photometry and res¬ 
olution, mid-IR studies intended to study the stellar content of 
galaxies beyond ~5 Mpc with current archival data from Spitzer 
will be severely limited. 

Ideal candidates are those that are isolated and meet photo¬ 
metric criteria in the IR that have been previously successfully 
used to identify massive stars. Because we have deep optical 
photometry (from the HST images) for our IR selected candi¬ 
dates, we can explore the differences and similarities in select¬ 
ing candidates using IR photometry versus selecting candidates 
using optical photometry. 

In pioneering the effort to hnd red supergiants in nearby 
galaxies, Massey| ( 1998 1 showed that while both the V - R and 
B - V photometric colors change with effective temperature, 
B - V also changes with surface gravity, and can therefore be 
used to uniquely identify RSGs from foreground. Milky Way red 
halo stars. A similar approach was implemented by 
(|2013|l with B - V and V - I colors. 


2005| and references therein), 41 from the LMC ( |Bonanos et'af 
2009 1, 19 from the SMC ( |Bonanos et al.|20T0| ), and 30 from the 


Dwarf Irregular galaxies (dirrs) in the Local Group ([Levesque & 


Massey|2012[ Britavskiy et al.|20T4| Britavskiy et al. in prep.). 


The box outlines the region described in Grammer & Humphreys 
(2013|l where Barmby et al. (|2006|l show that cluster candidates 


exist. Interestingly, this region also is predominantly the locus 
of the massive stars such as sgB[e] fromjBonanos et ^ (2009 


|2010| ) (blue dots), and the list of LB Vs, LBV candidates, Feii 
emission line stars, and other (blue) supergiants in M31 and M33 
from Humphreys et al. ( [2014 1 shown as green dots. Comparing 
our candidates with the locations of the spectroscopically classi¬ 
fied massive stars from the literature, objects in the upper right 
of the plot are most likely RSGs, while those in the lower left 
are their bluer massive star counterparts. The average reddening 
(Kim et al. 20121 is shown as an arrow, meaning some candidates 


Humphreys (20131 extracted a theoretical supergiant sequence 
from the stellar models of jBertelli et al.j ( |1994| ) to further illus¬ 
trate the usefulness of such a color-color diagram in selecting 
massive star candidates. Ten of the 112 candidates visually in¬ 
spected in the HST images were not recovered in one or more of 
the hlters needed for plotting on this color-color diagram. Figure 
shows the 102 (of 112) candidates listed in Table with B, 
V, and I magnitudes on a color-color diagram of B - V versus 
y - /. In total, there are 49 candidates with rank 2 or lower plot¬ 
ted in red in Figure|^and 53 candidates with a rank higher than 
2 plotted in black. Of the 49 highly ranked visually inspected 
candidates, 24 (49%) lie in the region of the color-color dia¬ 
gram similar to the RSG candidates of [Grammer & Humphrey sj 
( 2013| l. Within uncertainties, another four strong candidates also 
lie within this region. Based on the location of our candidates on 
the similar plot of B - V versus V - /, we expect that ~50% of 
the strong (49 red points in Figure candidates we selected 
via near- and mid-IR photometry are truly RSGs. To be cer¬ 
tain, we need follow-up spectroscopy to conhrm this hypothesis. 
For comparison, Massey ( 1998[ l applied a photometric cutoff in 
the B - V versus V - R colors, then performed spectroscopic 
follow-up observations in order to conhrm the objects as RSGs 
and found that for the faintest candidates, the success rate of 
the optical photometric criteria in selecting RSGs was 82%. To 
compare our expected success rate of recovering RSGs with that 
from optical surveys, we must note a few caveats. First, the op¬ 
tical photometry used by Massey et al. ( [2009| l, for example, is 
much more accurate than the photometry from the HST data. 
This is simply due to the much greater distance to the stars in 
M83 than to the studied stars in Local Group galaxies. The sec¬ 
ond caveat relates also to the accuracy of the photometry, this 
time from Spitzer. Selection criteria are only as good as the ac¬ 
curacy of the photometry, and as can be seen in Figures]^ and 
within uncertainties, some RSG candidates may move in and out 
of the region encompassed by our selection criteria. Because of 
these caveats, our predicted RSG success rate is lower than that 
from optical observations. 


outside the described regions are simply reddened RSGs in the 
upper part of the plot, or reddened sgB[e] stars (for example) 
in the lower central part of the plot. The strong candidates (red 
dots) in the lower right part of the plot are not very red 'm B -V 
while being quite red in V - /. [Grammer & Humphreys ( [2013 1 
put forth the explanation that objects in the bottom right of the 
color-color diagram may contain blended objects or compact H ii 
regions. Indeed, at the limit of the HST resolution, ultra compact 
H II regions, which are typically connected to small clusters of 
stars ( Conti et al.|2008| l, may be contaminants. In a study of the 
clusters in the central helds of the HST data for M83, [Chandar| 
([2010 1 show that reddening may cause the V - I values of 


et al. 


candidate LB Vs to move out of where they are expected in the 


same way as seen in our plot. [Whitmore & the WFC3 Science 
Oversight Committee ([201 1[) a^Kim et al.[^20l2|l bring up the 


possibility that these objects may be the chance superposition of 
a red and a blue star. Spectroscopic follow-up is the only way to 
definitively determine the nature of these objects as either red¬ 
dened LB Vs, a superposition of two different kinds of stars or 


simply clusters of stars. It should also be mentioned that Massey 
et al. ([2009 [l discussed that 5 of 19 spectroscopically conhrmed 


RSGs in M31 from Massey[ (ll998[l would not have met their 


photometric criteria. They stressed that RSGs may exist in other 
regions of the color-color diagram, but reiterated the point that 
their optically photometric criteria have a superb success rate of 
> 82% ( [Massey|1998[ l at selecting RSGs based on conhrmation 
from follow-up spectroscopy. Again, this is higher than our ex¬ 
pected ~50% success rate based on IR photometry. 

How does selecting evolved massive stars via IR photome¬ 
try compare with the [Massey[ ( [1998[ l optical selection method? 
Plotted in Figure hy color and galaxy are the spectroscopi¬ 
cally conhrmed RSGs shown i n grey in Figure [9| Also p lotted 
is the supergiant sequence from Grammer & Humphreys [ ( [2013[ l 
and [Bertelli et al. (1994 1 for reference. Because optical surveys 
were performed hrst and are thus more numerous and typically 
have better resolution, a very small number of conhrmed RSGs 


come from the IR-selected candidates in dirrs (Britavskiy et al. 
[2014[ Britavskiy et al. in prep.), which are shown as red triangles 
outlined in black. The spectroscopically conhrmed IR-selected 
candidates in the dirrs in the Local Group lie in exactly the same 
region as those spectroscopically conhrmed RSGs selected via 
optical photometry, and are represented by downward pointing 


6 








































































































































Williams et al.: Massive Stellar Content of M83 


red triangles that are not outlined in black. Clearly, there are too 
few RSGs selected via the IR to draw strong conclusions, but 
these first few attest to the usefulness of selecting RSGs from IR 
observations. Looking further at the plot, there are a few inter¬ 
esting possibilities that will be confirmed or rejected with more 
observations. For the confirmed RSGs in the SMC in Figure [T0| 
are the six RSGs with lower B-V representative of a metallicity 
effect? These RSGs are all K-type, and the trend of earlier aver¬ 
age spectral type of RSGs with decreasing metallicity of the host 
galaxy has been shown by Levesque & Mass^ ( |2012| l. More ob¬ 
servations of RSGs in metal-poor galaxies are needed to confirm 
or refute this possible trend. The Local Group dirr RSGs do not 
yet have deep enough photometry to indicate if this trend is real. 
Another question raised by this plot is why do the LMC RSGs, 
Milky Way RSGs, and M83 RSGs suffer from so much scatter 
on this plot? 

In Figure we plot RSGs grouped by four spectral types: 
early K, late K, early M, and late M. The K spectral type RSGs 
appear less reddened, and show less scatter. This follows from 
the intrinsic bluer colors of K-type RSGs compared to M-types. 
However, K-type RSGs also appear to have less reddening than 
their later-type companions. This is likely due to the fact they 
come from, on average, lower metallicity galaxies. One interest¬ 
ing possibility for this scatter for the later-type RSGs may be 


related to the behavior of HV 11423 in the SMC. Massey et al. 
( |2007| l tracked HV 11423 through changes in its spectral type 
from KOI to M5 I. The star changed by 2 mag in V, with changes 
in spectral type attributed to both the change in effective tem¬ 
perature, and the creation and dissipation of dust ([Massey et al. 


2007jl. HV 11423 is an exceptional case, but it would be interest¬ 


ing to investigate via long term observations whether the scatter 
of later spectral types in Figure [TT] is connected with a cycle of 
production and destruction of dust. 

What kind of RSGs do we expect to find in M83? The aver¬ 


age spectral type for RSGs in the Milky Way is M2 (Levesque 
& Massey|2012|l. The disk metallicity of M83 was shown to be 


1.9x Solar by Gazak et al. ( 2014| ). Thus, we may expect the aver¬ 
age spectral type of an M83 disk RSG to be M2 or later. In addi¬ 
tion, we may expect the typical M83 RSG to experience greater 
mass loss rates, as mass loss rates tend to increase with later 
spectral type. 


5. Conclusions 


We present point source catalogs for a 15' x 15' region centered 
on M83 based on Spitzer mid-IR and near-IR photometry. From 
these catalogs, we have selected candidate massive stars based 
on their mid- and near-IR photometric properties. These candi¬ 
dates have been culled by checking against catalogs of known 
objects, and also by inspection on high resolution HST images. 
The remaining 49 objects are strong candidates for evolved mas¬ 
sive stars, and await follow-up spectroscopy to further determine 
their nature and quantify the success rate of this technique in de¬ 
tecting evolved massive stars. 

We make the full catalog of point sources available to the 
community in order, for example, for future researchers to char¬ 
acterize the progenitors of core-collapse SNe. The methodology 
used in this paper may be readily applied to other nearby (< 5 
Mpc) galaxies to investigate their evolved massive stellar popu¬ 
lations. High quality archival Spitzer images exist for all nearby 
(< 10 Mpc) galaxies courtesy of the SINGS and LVL surveys 
( Kennicutt et al.|2003| Dale et al.|2009| ). Archival HST images 
also exist for many of the nearby galaxies. Analysis of follow¬ 
up spectroscopy of candidates discussed here for M83 will be 


needed to further test this methodology. However, studies start¬ 
ing merely with Spitzer photometry have already proved to be 
successful in discovering evolved massive stars ( jBritavskiy et al. 
|2014[ Britavskiy et al. in prep.). The validation of the useful¬ 
ness of this IR method of determining massive star candidates is 
important, especially because in many cases, deep optical pho¬ 
tometry simply does not exist. This method may also provide an 
alternative to the previously employed optical methods to inves¬ 
tigate the massive stellar content of the Local Universe. 
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Table 1: Infrared Photometry of 4528 Point Sources in the field of M83 


RA(J2000) 

(deg) 

Dec(J2000) 

(deg) 

J 

(mag) 

a-j 

(mag) 

K, 

(mag) 

(mag) 

[3.6] 

(mag) 

O' 3.6 

(mag) 

[4.5] 

(mag) 

O' 4.5 

(mag) 

[5.8] 

(mag) 

O's.s 

(mag) ... 

204.16129 

-29.85840 

16.40 

0.02 

16.01 

0.03 

15.91 

0.03 

16.03 

0.07 

14.90 

0.07 

204.16831 

-29.90027 

17.42 

0.03 

16.52 

0.03 

17.02 

0.04 

17.27 

0.11 

17.09 

0.28 

204.17050 

-29.81191 

16.64 

0.02 

16.27 

0.05 

16.47 

0.05 

16.60 

0.06 

16.21 

0.11 

204.17440 

-29.86086 

18.31 

0.05 

16.65 

0.06 

16.73 

0.06 

16.99 

0.08 

13.31 

0.01 

204.18935 

-29.91848 

16.28 

0.01 

15.86 

0.02 

15.91 

0.04 

16.07 

0.05 

15.01 

0.04 

204.20567 

-29.92462 

17.81 

0.04 

16.80 

0.05 

16.77 

0.05 

17.01 

0.07 

16.31 

0.25 

204.21006 

-29.91604 

18.40 

0.05 

16.71 

0.05 

16.85 

0.05 

17.08 

0.07 

13.79 

0.03 

204.21704 

-29.83906 

17.41 

0.03 

16.40 

0.04 

16.70 

0.05 

16.84 

0.05 

13.46 

0.17 

204.22253 

-29.82444 

17.43 

0.05 

15.97 

0.06 

15.55 

0.04 

15.66 

0.06 

15.92 

0.76 

204.23417 

-29.94974 

19.66 

0.15 

16.94 

0.10 

16.37 

0.03 

16.50 

0.06 

15.82 

0.11 


This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and 

content. 


Table 2: FourStar Photometry Compared to 2MASS Photometry 


2MASS ID 

■fpSF 

(mag) 

•f2MASS 

(mag) 

Difference 

(mag) 

(J - A'slpourStar 
(mag) 

(J - Ais)2MASS 
(mag) 

Difference 

(mag) 

13372193-2947529 

15.761(10) 

15.694(58) 

0.067 

0.889(23) 

0.888(141) 

0.001 

13371693-2947174 

16.105(8) 

16.067(88) 

0.038 

0.896(21) 

0.850(188) 

0.046 

13371216-2948184 

16.235(15) 

16.192(99) 

0.043 

1.027(37) 

0.923(195) 

0.104 

13364314-2947289 

15.781(12) 

15.872(77) 

-0.091 

0.981(46) 

0.947(145) 

0.034 

13364385-2949437 

15.427(11) 

15.448(63) 

-0.020 

0.767(42) 

0.634(130) 

0.133 

13363909-2949390 

15.696(7) 

15.729(77) 

-0.033 

0.284(28) 

0.242(226) 

0.042 

13363400-2949388 

16.011(12) 

16.000(82) 

0.011 

0.708(41) 

0.852(164) 

-0.144 

13363098-2950437 

15.155(10) 

15.194(48) 

-0.039 

0.850(36) 

0.806(91) 

0.044 

13363542-2953303 

16.219(10) 

16.366(130) 

-0.147 

0.993(24) 

0.950(250) 

0.043 

13370045-2957177 

15.670(7) 

15.712(65) 

-0.042 

0.487(25) 

0.463(178) 

0.024 

13370527-2956231 

15.962(10) 

15.939(80) 

0.023 

1.076(26) 

1.057(146) 

0.019 

13372091-2949532 

15.234(9) 

15.249(47) 

-0.015 

0.844(19) 

0.884(89) 

-0.040 


Table 3: Massive Evolved Star Candidates in M83 from HST, Near-IR, and Spitzer Photometry 


RA(J2000) 

(deg) 

Dec(J2000) 

(deg) 

U 

(mag) 

o-u 

(mag) 

B 

(mag) 

a-B 

(mag) 

V 

(mag) 

a-y 

(mag) 

I 

(mag) 

a-i 

(mag) 

J 

(mag) ... 

204.27068 

-29.78285 









18.07 ... 

204.27253 

-29.78501 










204.17782 

-29.78529 










204.33232 

-29.78565 










204.25152 

-29.80064 

26.67 

0.37 

26.56 

0.18 

24.08 

0.07 

20.78 

0.03 

18.74 ... 

204.24779 

-29.80129 

26.90 

0.42 

25.79 

0.13 

23.43 

0.05 

20.68 

0.03 

18.46 ... 

204.27008 

-29.80718 

25.25 

0.19 

24.72 

0.08 

22.35 

0.04 

20.21 

0.03 

18.28 ... 

204.31808 

-29.80785 










204.27208 

-29.80829 



29.69 

0.90 

26.68 

0.22 

22.04 

0.04 

18.94 ... 

204.23951 

-29.80921 

26.16 

0.30 

27.44 

0.27 

25.65 

0.13 

21.79 

0.04 

18.66 ... 


This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and 

content. 
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Fig. 1; Magnitude versus uncertainty for the Spitzer [3.6] and [4.5] point sources. The [4.5] values are in red and show both the 
higher uncertainties and brighter limiting magnitude limit for this bandpass. 
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[ 3 . 6 ]-[ 4 . 5 ] 


Fig. 2: Color-magnitude diagram (CMD) for M83 based on Spitzer IRAC Bands 1 and 2 photometry. Regions empirically known 
to contain red supergiants (RSGs) and supergiant Be stars (sgB[e]) in |Bonanos et al.|p009 |2010) l are outlined with red and blue 
boxes, respectively. Objects that lie within those regions are plotted with red and blue points. On the left side of the plot are shown 
the average Icr uncertainty ellipses for stars in 1 mag bins for the [3.6] photometry and the [3.6]-[4.5] color. For reference, also 
plotted are the locations of the prominent evolved massive stars rj Car ([Humphreys & Davidson 19941, M33 Variable A ([Humphrey 


et al. 2006 1 , and Object X (Khan et al. 20101 when placed at the distance (~4.8 Mpc) of M83. 
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Fig. 3; This figure is the same as Figurej^but includes grey points from the Spitzer Deep Wide Field Survey representing the expected 
regions of contamination from non-M83 sources. Note that only the more luminous and bluer RSG candidates lie distinctly outside 
the contaminated region. For reference, also plotted are the locations of the prominent evolved massive stars rj Car ( [Humphreys &| 
Davidson|1994i, M33 Variable A ([Humphreys et al.|2006[). Object X (Khan et al.|201 1 1 , when placed at the distance (~4.8 Mpc) of 


M83. 
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15 

14 

13 

12 

11 

10 



Fig. 6; Color-magnitude diagram of [3.6] vs. 7-[3.6] for targets in the Spitzer and near IR combined catalog. Regions of specific 
classes of stars as seen in Figure 3 of Bonanos et al. (2009 2010|l are outlined by a red box for RSGs, and a blue box for sgB[e] 


stars. 
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(e) Quality Rank 5 


Fig. 7: Composite HST images (from F336W U, F438W B, F547M y, F555W V, F657N Ha, F814W I) of objects from Table |3] 
demonstrating the “Quality” ranks listed in column 15. Both the F336W and F438W hlter images make up the blue colors, F547M 
or F555W images are shown as green, while red colors are from F657N Ha and F814W. Each image is 10" on a side. Over-plotted 
is the T.'7 Spitzer point spread function. The quality ranks are: (a) an isolated, excellent red candidate, rank 1, (b) an extremely 
red candidate, rank 2, (c) an unresolved cluster or other non-point source, rank 3, (d) a compact association, rank 4, (e) a likely 
background galaxy, rank 5. In these images. North is up and East is left. 
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204.35 204.30 204.25 204.20 204.15 


RA (J2000) (Degrees) 


Fig. 8; A 15' x 15' Spitzer 3.6 pm Band 1 image of M83, depicting objects with rank 1 and 2 as red circles, and rank 3 and higher 
shown as blue circles. All candidates lie within the HST coverage and are in the spiral arms, where massive star candidates are most 
likely to reside. 
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Fig. 9; Color-color diagram for M83 stars with HST photometry. Objects determined to be good candidates for follow-up spec¬ 
troscopy with rank 1 and 2 are shown in red, while objects with ranks 3 and higher are shown in black. Also plotted is a “typical” 
F-band reddening for stars in M83’s central region as reported by Kim et al. (2012|l. Spectroscopically confirmed RSGs in the Local 


Group are collectively shown as grey points. The dashed line represents the theoretical supergiant sequence from Bertelli et al. 
( 1994|l as used in|Grammer & Humphreys (|2013 i. The box outlines the region where cluster candidates are known to exist (Barmby 


|et al.|20()^ . The blue points are sgB[e] stars from the LMC and SMC ( |Bonanos et al.|200^|2010) l, while the green points are blue 
supergiants from M31 and M33 (Humphreys et al. 2014 1 . 
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Fig. 10; A color-color diagram showing the locations of spectroscopically confirmed RSGs broken down by galaxy in the Local 
Group. Over-plotted is the theoretical supergiant sequence from Grammer & Humphreys (20131 and Bertelli et al. (19941, as well 
as an arrow showing the typical reddening in M83 as determined by Kim et al. (2012|l. 
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Fig. 11; A color-color diagram showing the locations of spectroscopically confirmed RSGs in the Local Group broken down by 
spectral type. Over-plotted is the theoretical supergiant sequence from Grammer & Humphreys (2013|l and Bertelli et al. ( 1994|l. 
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(d) 204.31148-29.84119 

(e) 204.22728-29.86241 

(f) 204.33467-29.86917 



(g) 204.33385-29.87635 (h) 204.2828-29.89583 (i) 204.24671-29.90021 



G) 204.30452-29.90964 (k) 204.30305-29.91371 (1) 204.25023-29.91497 


Fig. 12: Objects from Tablej^with a rank of 1, thus excellent candidates for a massive star. The color scheme for these images is the 
same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 1'.'7 shown as a white circle. North is 
up and East is left. 
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(a) 204.21006-29.91601 (b) 204.22586-29.91788 (c) 204.21403-29.9242 




(g) 204.23951-29.80921 (h) 204.26044-29.81253 (i) 204.26763-29.81998 



(j) 204.26145-29.82107 (k) 204.25092-29.82219 (1) 204.2508-29.8277 


Fig. 13: More objects from Tablej^with a rank of 1, thus excellent candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 1"7 shown as a white circle. North 

is up and East is left. 
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(g) 204.25766-29.87868 (h) 204.26906-29.88573 (i) 204.25522-29.88789 



Fig. 14: More objects from Tablej^with a rank of 1, thus excellent candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 177 shown as a white circle. North 
is up and East is left. 
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(d) 204.24474-29.92718 


Fig. 15: More objects from Tablej^with a rank of 1, thus excellent candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 177 shown as a white circle. North 
is up and East is left. 
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(a) 204.24779-29.80129 (b) 204.27008-29.80718 (c) 204.27208-29.80829 



(d) 204.26209-29.81663 (e) 204.26922-29.81975 (f) 204.23924-29.84713 



Fig. 16: Objects from Tablej^with a rank of 2, thus decent candidates for a massive star. The color scheme for these images is the 
same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 1"7 shown as a white circle. North is 
up and East is left. 
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(d) 204.28307-29.84473 (e) 204.30836-29.84737 (f) 204.25821-29.85125 



(g) 204.31955-29.85543 (h) 204.31419-29.85995 (i) 204.23301-29.87065 



0)204.20346-29.87159 


Fig. 17: More objects from Tablej^with a rank of 2, thus decent candidates for a massive star. The color scheme for these images is 
the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of T.'7 shown as a white circle. North 
is up and East is left. 
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(a) 204.28665-29.81637 (b) 204.28041-29.81667 (c) 204.28956-29.81703 



(d) 204.25643-29.81915 (e) 204.26219-29.81984 (f) 204.26318-29.81985 



(g) 204.27432-29.82038 (h) 204.24071-29.83772 (i) 204.26684-29.84411 



G) 204.29662-29.85237 (k) 204.23938-29.86178 (1) 204.27568-29.87186 


Fig. 18: Objects from Tablej^with a rank of 3, thus not good candidates for a massive star. The color scheme for these images is the 
same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 1'.'7 shown as a white circle. North is 
up and East is left. 
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(a) 204.30318-29.87883 (b) 204.27873-29.88957 (c) 204.22592-29.89037 




(j) 204.29721-29.91526 (k) 204.29463-29.91608 (1) 204.27929-29.91728 


Fig. 19: More objects from Tablej^with a rank of 3, thus not good candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of T.'7 shown as a white circle. North 
is up and East is left. 
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(a) 204.26318-29.92358 (b) 204.26171-29.9246 (c) 204.24593-29.92507 



(d) 204.23161-29.89652 


Fig. 20: More objects from Tablej^with a rank of 3, thus not good candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of T.'7 shown as a white circle. North 
is up and East is left. 
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(a) 204.26152-29.81385 (b) 204.22893-29.8338 (c) 204.23292-29.83584 



(d) 204.24608-29.83645 (e) 204.24191-29.83685 (f) 204.30085-29.85674 



Fig. 21: Objects from Table with a rank of 4, thus very poor candidates for a massive star. The color scheme for these images is 
the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of T.'7 shown as a white circle. North 
is up and East is left. 
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(a) 204.2652-29.81649 (b) 204.23216-29.84833 (c) 204.21583-29.85924 



(0) 204.27851-29.8657 (e) 204.20375-29.89278 (1)204.25955-29.92432 

Fig. 22: More objects from Tablej^with a rank of 4, thus very poor candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of T.'7 shown as a white circle. North 
is up and East is left. 
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(a) 204.25765-29.81874 (b) 204.27278-29.86357 (c) 204.29479-29.90216 



(d) 204.22135-29.90849 


Fig. 23: Objects from Tablej^with a rank of 5, thus extremely poor candidates for a massive star. The color scheme for these images 
is the same as for Figure]^ Each image is 10" on a side, with the Spitzer point spread function of 1'.'7 shown as a white circle. North 
is up and East is left. 
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